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For  the  numerical  simulation  of  Indian  Ocean  monsoon  de¬ 
pression,  an  idealized  zonal  current  pertinent  to  summer 
monsoon  conditions  is  used  as  the  basic  flow.  Initial  perturba¬ 
tions  superposed  on  the  basic  flow  consist  of  a  prescribed 
heating  source.  The  simulated  winds  at  4S  and  72  h  are  dis¬ 
cussed.  The  simulation  reproduces  many  observed  phenomena, 
including  the  storm  track.  But  the  cumulus  parameterization 
currently  used  in  the  NRL  mesoscale  model  yields  unrealistically 
heavy  rainfall  on  steep  mountain  slopes.  It  needs  modification 
in  future  use. 

Numerical  simulation  is  also  carried  out  to  simulate  the 
interaction  of  two  tropical  cyclones  (Fujiwhara  effects). 

Results  show  that  on  an  f -plane  with  no  mean  wind,  the  move¬ 
ments  of  the  two  cyclones  consist  of  a  mutual  cyclonic  rotation, 
mutual  attraction,  and  eventual  coalescense.  The  trajectories 
of  two  tropical  cyclones  of  equal  (different)  strength  are 
symmetric  (asymmetric)  relative  to  the  center  of  mass  of  the 
two  storms.  Further  numerical  experiments  indicate  that  the 
latitudinal  variation  of  the  Coriolis  parameter  has  great 
influence  on  the  trajectories. 

Concerning  easterly  waves,  we  attempted  to  simulate 
the  transformation  of  an  easterly  wave  into  an  incipient 
tropical  depression.  In  a  uniform  easterly  current,  the  NRL 
mesoscale  model  successfully  simulates  the  formation  of  an 
easterly  wave.  It  has  very  realistic  latitudinal  components 
and  relative  humidity  field,  but  the  upward  motion  in  the 
trough  is  too  strong.  We  suggest  for  future  study  latitudinal 
shear  in  the  mean  current  and  variation  of  the  Coriolis  parameter 
be  included  in  the  model. 
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!  I  CHAPTER  1.  INTRODUCTION 

I 

I 

This  final  technical  report  documents  the  research  efforts 
I  in  the  numerical  simulation  of  tropical  weather  systems  carried  out 

i 

under  the  auspices  of  ONR  Contract  N00014-81-C-2224,  for  the  Atmos¬ 
pheric  Physics  Branch,  Environmental  Sciences  Division,  of  the 
,  Naval  Research  Laboratory  (NRL) . 

The  research  efforts  can  be  topically  separated  into  three 
'  areas;  namely,  monsoon  depressions,  interactions  of  typhoon  pairs, 

and  easterly  waves.  In  the  following  three  chapters,  we  will  dis¬ 
cuss  In  detail  the  research  efforts  and  results  obtained  in  each 
of  these  areas.  In  short,  the  following  tasks  have  been  completed: 

(1)  A  dynamic  initialization  scheme  through  prescribed  heating  was 
devised  so  that  dynamically  consistent  tropical  disturbances  Including 
the  convergence  component  can  be  realistically  simulated  in  the  NRL 

♦ 

mesoscale  model. 

(2)  Major  modifications  of  the  NRL  mesoscale  model  were  carried 
out  to  facilitate  simulation  studies  of  tropical  weather  systems. 

> 

The  modifications  include  changes  in  grid  system  structure,  spatial 
resolution,  cyclic  boundary  conditions,  and  incorporation  of  topo¬ 
graphy. 

i 

(3)  Realistic  simulations  of  easterly  waves  and  monsoon  depressions 
have  been  accomplished.  Problem  areas  for  future  research  have  been 
identified  and  solutions  are  suggested. 


1 


In  addition,  all  computer  codes  generated  in  connection  with  our 
efforts  have  been  transferred  to  Naval  Research  Laboratory  personnel 
and  are  stored  in  the  Tl-ASC  in  both  softcopy  and  hardcopy  forms.  Some 
of  these  codes  can  be  implemented  in  the  Navy's  operational  forecasting 
models  to  improve  weather  forecasting  in  a  number  of  tropical  regions. 


I 


CHAPTER  2.  MONSOON  DEPRESSION 


2.1  BACKGROUND 


The  summer  and  winter  monsoons  of  south  Asia  constitute  very  pro¬ 
nounced  regional  anomalies  in  the  general  circulation  of  the  atmosphere 
resulting  from  land-sea  contrasts  and  orographic  features.  The  land 
masses  of  Africa  and  Eurasia  cover  about  60%  of  the  northeast  quadrant 
of  the  globe  whereas  nearly  80%  of  the  area  is  covered  by  water  in  the 
other  three  quadrants.  Unlike  the  Atlantic  and  Pacific  Oceans,  the  In¬ 
dian  Ocean  is  landlocked  to  the  north  by  the  Asian  continent.  The  geo¬ 
graphic  setting  gives  rise  to  extreme  thermal  contrast  between  the  land 
and  the  sea  in  the  summer  and  in  the  winter,  which  is  the  crucial  fac¬ 
tor  in  the  development  of  monsoon  circulations  in  this  region. 


The  mountain  systems  in  this  region  play  an  Important  role  in 
the  generation  and  maintenance  of  the  south  Asian  summer  monsoon.  The 
important  features  include  the  Himalayas  and  Tibetan  plateau  to  the 
north,  the  Western  Ghats  on  the  Indian  Peninsula,  and  the  relatively 
flat  Gangetic  valley.  The  sharp  gradient  at  the  southern  edge  of  the 
Tibetan  plateau  inherently  makes  numerical  simulation  for  this  region 
difficult. 
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Our  study  is  concerned  with  the  numerical  simulation  of  monsoon 
depressions  that  frequently  occur  during  the  summer  monsoon  season. 

Monsoon  depressions  usually  form  over  the  warm  Bay  of  Bengal  and 
move  morthwestward  inland,  bringing  heavy  precipitation  to  the  northern 
Indian  peninsula.  These  depressions  can  produce  strong  winds  and  weather 
phenomena  in  the  northern  Indian  Ocean  and  the  northern  Arabian  sea. 

It  is  hoped  that  sucessful  numerical  simulations  of  these  depressions  can 
lead  to  a  better  understanding  of  the  relevant  physical  processes 
causing  their  genesis,  intensification,  decay,  and  interactions  with 
larger  scale  circulation  in  this  region.  Such  knowledge  will  eventually 
lead  to  better  regional  weather  forecasts. 

2.2  MODIFICATION  OF  THE  MESOSCALE  MODEL 

Due  to  the  high  mountains  and  their  proximity  to  the  model  lateral 
boundary  it  is  nearly  impossible  to  maintain  a  zero  net  mass  inflow, 
which  is  one  of  the  essential  requirements  of  the  model,  within  the  domain 
shown  in  Figure  2.1.  Thus  ten  grid  points  were  added  to  each  of  the 
east  and  west  lateral  boundaries.  The  surface  elevations  of  these  grid 
points  are  linearly  interpolated  in  the  x-direction  so  that  the  elevations 
are  uniformly  zero  at  the  outermost  grid  points.  Gradually  increased 
kinematic  viscosity  coefficients  are  applied  to  these  additional  grid 
points  to  serve  as  sponge  boundaries.  Because  the  prevailing  large 
scale  winds  are  predominately  easterlies  or  westerlies,  no  additional 
grid  points  are  needed  along  the  north  and  south  boundaries.  However, 
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increased  viscosity  coefficients  are  applied  there  too.  There  are  a  total 
of  71  x  31  horizontal  grid  points  and  five  sigma  layers  in  the  vertical. 


with  1°  horizontal  resolution. 

2.3  PRELIMINARY  RESULTS 

An  idealized  zonally  mean  field  similar  to  that  observed  is  used 
as  the  basic  flow  field  (Fig.  2.2).  The  mean  zonal  flow  includes 
westerlies  at  low  levels  with  maximum  speed  of  V7  m  s-1  at  13°  and 
easterlies  at  high  levels  with  a  jet  of  more  than  25  m  s-1  at  low  lat¬ 
itudes  (Krishnamurtl,  1971).  In  geostrophic  balance  with  the  mean  flow, 
the  temperature  field  has  a  10  K  north-south  gradient  in  the  midtroposphere. 

The  perturbations  which  initiate  the  monsoon  depression  are  intro¬ 
duced  to  the  model  as  prescribed  heating  pertinent  to  a  cloud  c-uster. 

Saha  et  al  (1981)  found  that  more  than  75%  of  monsoon  depressions  that 
occur  in  the  Bay  of  Bengal  can  be  traced  by  surface  pressure  tendency 
back  to  typhoons  or  tropical  depressions  in  the  Pacific  Ocean.  The 
remnants  of  these  tropical  disturbances,  posing  as  cloud  clusters, 
can  redevelop  once  they  move  over  the  very  warm  Bay  of  Bengal.  Earlier 
attempts  were  made  to  initiate  monsoon  depressions  in  our  model  by  pre¬ 
scribing  pressure  tendency  patterns.  These  attempts  failed  to  generate 
strong  enough  perturbations  in  the  momentum  field,  obviously  due  to  the 
unfavorable  mass-momentum  adjustments  on  the  scale  of  monsoon  depressions. 
The  prescribed  heating  currently  applied  produces  suitable  response  in 
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the  momentum  field.  The  prescribed  heating  in  the  preliminary  test 
presented  next  is  kept  stationary  for  24  h.  The  heating  in  the  model 
after  24  h  is  determined  by  the  cumulus  parameterization  scheme. 

Fig  2.3  shows  the  surface  (c  *  0.925)  flow  pattern  at  48  h,  when 
the  center  of  the  depression  has  just  moved  inland.  The  initial  west¬ 
erly  jet  at  M3°  is  enhanced  over  the  Indian  peninsula,  and  is  turning 
into  a  southwesterly.  It  forms  a  trough  in  the  Bay  of  Bengal,  which  sup¬ 
plies  ample  water  vapor  into  the  depression.  There  is  a  directional  wind 
shear  north  of  the  depression  caused  by  the  gradient  of  the  topography 
on  the  southern  edge  of  the  Tibetan  plateau.  The  maximum  wind  shear 
near  the  depression  center  is  ^25  m  s-*,  in  agreement  with  the  obser¬ 
ved  strength.  The  depression  moves  west-northwest,  against  the  surface 
westerlies,  also  in  agreement  with  observations  (Sikka,  1978). 

It  appears  there  are  two  major  air  streams  into  the  depression.  One 
extends  from  the  westerly  jet,  and  the  other  over  the  western  Gangetic 
valley  extends  northwestward  from  the  depression  center.  This  suggests 
that  the  depression  has  frontal  structure .  Besides  the  depression  center, 
there  are  two  areas  with  heavy  precipitation.  One  area  is  along  the  wind- 
slopes  along  the  west  coast  of  the  Indian  Peninsula  between  7-17°  lat¬ 
itude.  The  heavy  and  concentrated  precipitation  there  is  due  to  the  lift¬ 
ing  and  frictional  effects  of  the  western  Ghats.  Heavy  rainfall  during 
the  summer  monsoon  is  often  observed  there  (Ananthakrishnan,  1978) .  The 
second  area  of  heavy  precipitation  outside  the  depression  is  associated 
with  a  strong  convergence  zone  along  the  steep  mountain  slope  at  the 
southeast  edge  of  the  Tibetan  plateau  near  80°E  and  27°N.  In  fact. 
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this  "convergence"  is  fictitious.  It  arises  because  the  lowest  sigma 
layer  follows  the  terrain  in  our  model,  in  which  vertical  wind  shears 
appear  as  horizontal  convergence.  This  fictitious  convergence  causes 
an  erroneous  release  of  latent  heat  through  the  cumulus  parameterization. 
Part  of  our  future  effort  will  focus  on  the  improvement  of  Kuo's  scheme 
in  its  application  to  sloped  terrains. 

Fig  2.4  shows  the  surface  wind  field  at  72  h.  The  surface  center 
has  moved  inland  in  a  west-northwest  direction  with  the  maximum  wind 
slightly  decreased  to  22  m  s-*.  The  decreased  intensity  is  mainly  due 
to  increased  surface  friction  and  lack  of  moisture  once  the  fetch  from 
warm  ocean  becomes  far.  The  Tibetan  plateau  has  a  great  influence  on 
the  southly  air  flow  originated  from  the  Bay  of  Bengal  and  there  appears 
to  be  a  stagnant  point  at  90°E  and  26°N.  The  surface  winds  are  parallel 
to  the  slope  west  of  this  point. 


2.4  SUMMARY 


The  initiation  and  movement  of  an  idealized  monsoon  depression 
has  been  simulated  with  the  NRL/SAI  numerical  model.  The  model  domain 
covers  an  area  from  50°  -  100°E  and  5°  -  35°N.  To  improve  the  model, 
ten  grid  points  have  been  added  to  the  east  and  west  lateral  boundaries 
to  smooth  out  the  terrain  and  to  serve  as  a  sponge.  There  are  five  layers 
in  the  vertical,  and  the  horizontal  resolution  is  one  degree. 
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An  idealized  zonal  current  pertinent  to  summer  monsoon 
conditions  is  used  as  the  basic  current.  Perturbations  superposed 
on  the  mean  current  consist  of  a  prescribed  heating  source  appro¬ 
priate  to  a  propagating  cloud  cluster.  Simulated  surface  wind 
fields  at  48  and  72  h  are  presented  and  discussed.  The  simulation 
reproduces- many  observed  phenomena  associated  with  monsoon  depressions. 
But  Kuo's  cumulus  parameterization  is  found  to  give  unrealistic 
heating  near  mountain  slopes  where  vertical  wind  shear  exist. 

Its  application  in  the  model  in  such  regions  needs  further  mod¬ 
ification. 
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CHAPTER  3.  THE  INTERACTIONS  BETWEEN  TWO  TROPICAL  CYCLONES 


3.1  BACKGROUND 


When  two  tropical  cyclones  are  present  simultaneously  in 
the  same  region,  it  is  often  observed  that  they  move  around 
each  other  cyclonically.  The  phenomenon  was  first  mentioned 
in  Fujiwhara  (1921),  in  which  a  mathematical  approach  was  un¬ 
dertaken  where  the  typhoons  were  represented  by  solid  rotating 
vortices.  Haurwitz  (1951),  examined  several  tropical  cyclone 
pairs  by  introducing  the  concept  of  the  center  of  mass  around 
which  the  two  tropical  cyclones  rotate  at  a  rate  given  by  the 
gradient  winds.  Liu  and  Wang  (1966),  studied  some  recent  cases 
in  the  Pacific  Ocean.  They  found  that  the  two  interacting  tropical 
cyclones  are  not  always  attracted  to  each  other,  especially 
when  there  is  a  large  shear  in  the  basic  flow. 

Because  most  tropical  cyclone  pairs  occur  over  the  open 
ocean  where  data  is  sparse,  researchers  had  to  work  with  no 
more  than  conventional  isobaric  analyses  with  no  details  of 
flow  fields.  These  analyses  may  contain  factors  such  as  shears 
in  the  basic  flow  and  changes  in  the  Intensity  that  can  obscure 
the  true  Interactions.  Our  effort  was  to  simulate  this  phenomenon 
with  a  numerical  model  and  to  isolate  the  Fu.liwhara  effects 
in  an  environment  with  no  large  scale  flow. 
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3.2  THE  NUMERICAL  MODEL 

The  numerical  model  is  identical  to  the  one  in  Chang  and  Madala  (1980), 
except  for  the  parameterization  of  the  latent  heating.  Here,  a  pre¬ 
scribed  heating  is  used.  The  heating  rate  is  defined  as 


i  &  7rr 

I  Q  COS  ( —  )  cos 
O  R 


Q  (r,o) 
o 


v (o— 0 . 5) 


0.45 


,  for  r£R  ,  O.HcK  0.9 


(3-1) 


0,  otherwise 

where  r  is  the  distance  between  a  grid  point  and  the  low  pressure  cen¬ 
ter,  R  «  300  km  is  the  radius  of  the  heating  pattern,  and  0  is  the  max- 

o 

mum  heating  rate.  Two  values  of  Q  ,  100  K  day-1  and  200  K  day-1,  are 

o 

used  in  various  experiments  to  define  the  weak  and  strong  tropical  cyc¬ 
lone.  The  spatial  distributions  of  the  heating  pattern  are  illustrated 
in  Fig.  3.1.  The  heating  pattern  follows  the  low  pressure  centers  which 
are  determined  by  a  linear  interpolation  of  surface  pressures. 

The  model  has  51  x  51  horizontal  grid  points  with  seven  sigma  lay¬ 
ers  in  the  vertical.  The  horizontal  resolution  is  1/2°.  The  east-west 
boundaries  are  cyclic,  and  a  sponge  boundary  condition  is  applied  to  the 
north  and  the  south  boundaries. 

The  heating  prescribed  by  (3.1)  generates  a  very  realistic  flow  field 

for  tropical  cyclones.  Figure  3.2  shows  the  radial  distribution  of  the 

quasi-steady  wind  speeds  at  the  6th  and  the  7th  model  layers  after 

24  h  of  stationary  heating  with  $  ■  200  K  day  -1 .  The  wind  speeds  have 

o 

a  peak  at  r  -  1°  and  decrease  gradually  outward  without  any  noticeable 
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discontinuity  at  r»R.  We  note  that  by  using  the  prescribed 
heating,  we  cannot  simulate  the  effects  of  the  interaction  between 
the  two  storms  on  the  cumulus  convection  in  each  individual 
storm,  which  in  turn  may  alter  the  interactions  by  changing 
the  storm's  intensity.  But  these  effects  may  be  secondary  and 
only  important  when  the  tropical  cyclone  pair  are  close.  As 
a  preliminary  study,  we  opt  to  use  the  more  economical  prescribed 
heating  to  study  the  first  order  effects. 

The  two  tropical  cyclones  are  dynamically  initialized  by 
24  h  of  stationary  heating  at  two  locations  M0°  in  longitude 
apart  and  at  the  same  latitude.  Fig  3.3  shows  the  surface  pressure 
field  at  24  h  in  a  case  where  strong  heating  (Qq  ■  200  K  day"1) 
is  applied  at  both  locations.  The  southward  displacement  of 
the  western  storm  (storm  A)  and  the  northward  displacement  of 
the  eastern  storm  (storm  B)  indicate  that  their  interaction 
has  already  caused  the  two  storms  to  rotate  cyclonically ,  in 
spite  of  the  stationary  heating.  Note  the  heating  is  allowed 
to  move,  following  the  low  pressure  centers  after  the  24  h  ini¬ 
tialization  period.  All  numerical  experiments  in  our  study 
are  listed  in  Table  3.1. 


TABLE  3.1 


» 

EXP. 

0  for  storm  A 

0 

(K  day-1) 

6  for  storm  B 

0 

(K  day-1) 

Coriolis 

Parameter  (s 

1 

1. 

200 

200 

f«4.37xl0-5 

2. 

100 

100 

»r 

3. 

100 

200 

it 

1 

4. 

200 

200 

real  f 

5. 

100 

200 

fl 

6. 

200 

100 

It 

1 

7. 

200 

“ 

M 

3.3.1  ON  f-PLANE 

Exps.  1-3  are  Intergrated  with  a  constant  f.  Therefore  the  orientation 
has  no  meaning  and  the  results  are  independant  of  the  absolute  initial  pos¬ 
itions  of  the  storms.  Figure  3.4  shows  the  trajectories  of  the  storm  cen¬ 
ters  in  Exp.  1,  in  which  two  model  tropical  cyclones  are  of  the  same 
strength.  The  trajectories  show  that  the  two  storms  rotate  about  each 
other  in  a  cyclonic  fashion  before  the  coalescence  after  96  h.  The  two 
trajectories  are  symmetric  about  a  center  of  mass,  which  happens  to  be  the 
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center  of  the  model  domain.  Superimposed  on  the  symmetric  rotation 


there  is  a  convergence  of  the  two  tropical  cyclones.  The  distance 
between  the  two  storms  decrease  from  M024  km  at  24  h  to  ^612 
km  at  96  h.  The  symmetry  remains  until  after  96  h  when  the 
two  heating  patterns  overlap  and  one  single  large  area  of  low 
pressure  is  formed. 

Exp.  2  is  identical  to  Exp.  1  except  that  the  heating  rate 
is  halved.  The  cyclonic  trajectories  (Fig.  3.5)  are  still  re¬ 
markably  symmetric  about  the  center  of  mass.  Because  of  the 
weaker  heating,  two  identifiable  centers  still  exist  at  120  h 
when  they  are  only  MOO  km  apart.  We  take  note  again  that  impor¬ 
tant  feedbacks  between  cumulus  convection  and  the  storm  pair's 
interaction  may  have  been  lost  due  to  the  prescribed  heating 
in  our  model. 

The  speeds  at  which  the  two  tropical  cyclones  in  Exps, 

1  and  2  rotate  around  and  approach  each  other  are  shown  in  Fig. 
3.6.  The  tangential  velocity  of  the  cyclonic  rotation  in  Exp. 

1  increases  from  M  m  s-1  at  24-36  h  to  well  over  6ms-1  after 
72  h  as  the  separation  between  the  two  tropical  cyclones  becomes 
small.  The  rotation  speeds  in  Exp.  2  are  about  1  m  s-1  slower 
than  those  in  Exp.  1.  Haurwitz  (1951),  pointed  out  that  the 
rotation  rate  is  proportional  to  the  combined  circulation  strength 
of  the  two  tropical  cyclones,  which  was  defined  as  the  mass 
flow  between  two  subjectively  selected  Isobars,  Using  the  average 
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maximum  winds  in  our  model  as  a  measure  of  the  combined  strength, 
we  found  the  relationship  between  the  rotation  rate  and  the 
strength  of  the  interacting  tropical  cyclones  postulated  by 
Haurwitz  (1951)  to  be  neither  clear  nor  well-defined.  It  is 
surprising,  however,  that  the  rate  of  convergence  is  independent 
of  the  combined  strength  as  indicated  by  the  radial  velocities 
in  Fig.  3.6. 

Fig.  3.7  shows  the  trajectories  of  the  two  storm  centers 
in  Exp.  3,  in  which  the  maximum  heating  rate  in  storm  A  is  only 
half  that  of  storm  B.  We  see  that  the  two  storms  still  rotate 
about  each  other  cyclonically  and  that  they  still  move  toward 
each  other.  However,  the  trajectories  are  asymmetric  and  the 
weaker  storm  A  moves  much  faster  than  the  stronger  storm  B  in 
a  way  similar  to  that  of  a  binary  celestial  system  in  which 
the  two  bodies  have  different  masses.  This  type  of  interaction 
has  been  observed  between  Typhoons  Flossie  and  Grace  (1950) 

(Liu  and  Wang  1966).  The  center  of  rotation  is  not  stationary. 
Instead,  it  moves  within  a  small  area  defined  by  the  lines  con¬ 
necting  the  two  storm  centers  at  different  times. 

Results  in  Exps.  1-3  demonstrate  that  on  an  f-plane  the 
interactions  between  two  tropical  cyclones  consist  of  a  cyclonic 
rotation,  attraction,  and  eventual  coalescence. 
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In  order  to  examine  the  airflow  fields  associated  with 
the  interaction,  we  transform  the  model  500  mb  wind  fields  in 
Exp.  1  onto  a  polar  grid  with  respect  to  the  center  of  the  model 
domain.  We  now  define  the  azimuthal  mean  velocity  as 


wher-s  v  ,  v^,  are  radial,  tangential  and  vertical  velocities 
on  the  polar  grid,  and  0  is  the  azimuthal  angle.  Fig.  3.8  shows 
the  azimuthally  mean  vertical,  tangential,  and  radial  velocities 
at  24,  48,  72,  and  96  h.  It  is  Interesting  that  the  mean  momentum 
fields  relative  to  the  center  of  domain  shown  here  are  very 
similar  to  those  in  intensifying  tropical  cyclones.  For  example, 
at  24  h  there  is  a  maximum  mean  tangential  velocity  of  ^  m  s-1 
at  r  ^600  km  and  a  minimum  of  ^-2  ms-1  at  r  ^400  km  reflecting 
the  cyclonic  wind  fields  about  the  two  storm  centers.  The  maximum 
tangential  velocity  gradually  increases  to  M6  m  s”*  and  moves 
toward  the  center  to  ^>350  km  at  96  h.  The  maximum  inflow  also 
develops  from  2ms-1  at  r  n* 600  km  at  24  h  to  'U  m  s"1  at  r  ^350  km 
at  96  h.  The  evolution  of  the  mean  vertical  velocity  is  also 
similar  to  that  of  an  intensifying  and  contracting  tropical 
cyclone. 

The  development  of  the  azimuthal  mean  circulation  can 
also  be  illustrated  by  comparing  the  kinetic  energy  (KE)  of 
the  mean  velocity  (KEM)  and  the  KE  of  the  eddy  velocity  (KEE) , 

Where 
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KEM  • ffu 2  [v  \]  •  [v  i  r  dr  d«  , 

KEE  ~JJ 1/2(V,  v't)  .  (v^,  ♦’)  r  dr  de  ,  (3.3) 

where  (v^,  v*)  -  ( v_,  vt)  - 

As  shown  by  Fig.  3.9,  the  KEE,  which  can  mostly  be  attributed 
to  the  circulations  around  the  two  centers,  reaches  a  quasi¬ 
steady  state  after  36  h.  Meanwhile  the  KEM,  representing  the 
strength  of  the  mean  circulation  as  depicted  by  Fig.  3.8  around 
the  center  of  rotation,  steadily  increases  until  the  coalescence 
of  the  two  tropical  cyclones.  The  ratio  of  KEE/KEM  decreases 
from  ^3  at  24  h  to  less  than  1  at  96  h. 

These  analyses  suggest  that  a  mean  circulation  relative 
to  the  center  of  rotation  develops  due  to  the  interaction  of 
two  tropical  cyclones.  This  mean  circulation  includes  tangential 
and  radial  components  resembling  those  associated  with  tropical 
cyclones.  The  cyclonic  and  convergent  movements  of  two  inter¬ 
acting  tropical  cyclones  are  believed  to  be  caused  by  the  ad- 
vective  effects  of  the  "mean"  circulation. 

3.3.2  THE  EFFECTS  OF  VARIATIONS  OF  THE  CORIOLIS  PARAMETER 

Exps.  4-7  were  carried  out  with  variable  Coriolis  parameter. 
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which  can  produce  a  northwestward  drift  of  tropical  cyclone 
in  the  northern  hemisphere  (Adem,  1956;  Madala  and  Piascek, 

1975;  Chang  and  Madala  1982;  Holland,  1982).  The  velocity 
of  the  drift  depends  on  the  latitude  and  on  the  cyclone's  strength. 
To  examine  the  free  drift  of  a  single  tropical  cyclone  in  our 
model,  we  carried  out  Exp.  7.  As  shown  by  trajectory  C  in 
Fig.  3.10,  the  model  tropical  cyclone  has  an  initial  northward 
movement,  but  changes  toward  the  northwest  after  36  h.  The 
0-72  h  mean  drift  velocity  is  1.18  m  s“*  toward  the  west  and 
1.37  m  s"1  toward  the  north.  The  center  at  72  h  is  '^6°  to  the 
west  and  ^6.5°  to  the  north  of  the  initial  position. 

The  latitudinal  variations  of  the  Coriolis  parameter  have 
a  very  pronounced  effect  on  the  trajectories  of  the  two  inter¬ 
acting  tropical  cyclones.  The  trajectories  of  the  two  tropical 
cyclones  with  equal  strength  in  Exp.  4  is  shown  in  Fig.  3.10. 
Instead  of  rotating  around  the  +  point  as  in  Exp.  1,  the  tropical 
cyclone  to  the  west  (storm  A)  moves  toward  the  southeast  then 
quickly  turns  toward  the  northeast,  while  the  tropical  cyclone 
to  the  east  (storm  B)  rapidly  moves  northwestward  and  rotates 
cyclonically  with  respect  to  storm  A.  The  two  storms  eventually 
coalesce  at  87  h,  with  storm  B  having  traveled  a  much  larger 
distance  from  its  initial  position  than  storm  A.  We  have 
computed  the  trajectories  A  and  B  relative  to  storm  C(Exp. 

7)  in  an  attempt  to  remove  the  beta-drift  —  -although  the  removal 
is  Incomplete  in  that  storm  C,  while  having  the  same  strength, 
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cannot  be  at  the  same  latitude  and  therefore  the  same  drift 
speed  simultaneously  with  both  storms  A  and  B.  The  resultant 
relative  trajectories  (not  shown)  are  very  similar  to  those 
in  Exp.  1  (Fig.  3.3).  We  therefore  conclude  that  the  asymmetric 
trajectories  exhibited  by  two  interacting  tropical  cyclones 
on  a  real  f  surface  consist  of  the  symmetric  rotation  found 
on  an  f-plane  and  a  beta  drift  toward  the  northwest. 

Exp.  5  is  to  compared  with  Exp.  3,  where  storm  A  is  weaker 
than  storm  B.  The  trajectories  of  the  storm  centers  in  Exp. 

5  (Fig.  3.11)  again  appear  very  different  from  those  in  Exp. 

3.  The  stronger  storm  B  shows  more  notlcable  northwest  drift 
than  in  Exp.  3.  The  weaker  storm  A  rotate's  cyclonically  toward 
the  southeast  at  a  much  reduced  rate  and  with  a  smaller  radius, 
apparently  due  to  the  counteracting  beta  drift. 

Most  interesting  is  Exp.  6,  in  which  storm  A  is  stronger 
than  storm  B.  From  24  to  72  h,  the  weaker  storm  B  moves  cycloni¬ 
cally  relative  to  storm  A.  In  the  meantime  storm  A  moves  slowly 
toward  the  southeast  nearly  perpendicular  to  and  away  from 
storm  B.  We  believe  that  this  is  due  to  a  combination  of  the 
Interaction  and  the  stronger  beta  drift  of  storm  A.  The  trajec¬ 
tories  take  a  strange  turn  after  72  h  because  the  two  storms 
start  to  influence  each  other  through  the  cyclic  east-west 


boundaries 


It  can  be  concluded  that  the  latitudinal  variation  of 
the  Coriolis  parameter  has  profound  impact  on  the  trajectories 
of  two  interacting  storms.  The  northwest  drift  caused  by  the 
beta  effect  greatly  distorts  the  symmetric,  cyclonic  rotation 
found  on  an  f-plane.  The  differential  beta  drift  also  can 
make  two  interacting  tropical  cyclones  move  apart  if  the  storm 
to  the  west  is  stronger. 

3.4  SUMMARY 

The  Interaction  between  two  tropical  cyclones  (the  Fujiwhara 
effect)  is  simulated  and  studied  with  the  NRL/SAI  mesoscale 
numerical  model.  Included  in  the  model  are  parameterizations 
of  boundary  layer  and  latent  heating  with  presribed  spatial 
distribution  moving  with  the  low  pressure  center.  The  model 
contains  51  x  51  x  7  grid  points  with  1/2  degree  horizontal 
resolution.  The  east-west  boundary  conditions  are  cyclic  and 
the  north-south  boundaries  are  sponge.  The  two  tropical  cyclones 
are  dynamically  initialized  at  the  same  latitude  by  a  24  h 
stationary  heating  at  two  centers  10°  in  longitude  apart. 

Results  show  that  on  an  f-plane  with  no  mean  wind,  the 
movements  of  the  two  cyclones  consist  of  8  mutual  cyclonic 
rotation  and  attraction.  The  spiral-like  paths  of  the  two  cyclones 
of  equal  (different)  strength  are  symmetric  (asymmetric)  relative 
to  the  "center  of  mass"  of  the  two  cyclones.  The  rate  of  rotation 


is  approximately  proportional  to  the  combined  strength  of  the 
two  tropical  cyclones  while  the  rate  of  convergence  is  not. 

The  two  tropical  cyclones  eventually  coalesce. 

Further  analyses  show  that  the  cyclonic  rotation  and  the 
mutual  attraction  of  two  interacting  cyclones  are  related  to 
the  development  of  a  "mean"  circulation  about  the  center  of 
mass,  similar  in  itself  to  a  weak  tropical  cyclone,  while  the 
"eddy"  kinetic  energy  about  the  centers  of  the  two  interacting 
tropical  cyclones  remains  relatively  invariant. 

When  the  latitudinal  variation  of  the  Coriolis  parameter 
is  included  in  the  model,  the  movements  of  two  interacting 
tropical  cyclones  are  complex.  The  trajectories  of  two  equally 
strong  tropical  cyclones  include  a  symmetric  cyclonic  rotation 
with  a  northwestward  beta  drift.  The  rate  of  their  mutual 
rotation  seems  to  be  unaffected  by  the  variable  Coriolis  parameter. 
The  location  of  the  eventual  coalescence  is  ,v60  to  the  north¬ 
west  of  the  initial  center  of  mass.  If  the  tropical  storm 
to  the  west  is  stronger  than  the  one  to  the  east,  we  found 
that  they  move  away  from  each  other  because  the  differential 
beta  drifts  overcome  the  mutual  attraction. 
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CHAPTER  4  EASTERLY  WAVES 


4 . 1  BACKGROUND 


The  existence  of  westward  traveling  wave  disturbances 
in  the  easterly  tropical  belt  has  been  established  since  the 
studies  by  Riehl  (1945)  and  Palmer  (1952).  Later  studies  yielded 
a  more  dynamically  and  thermodynamically  consistent  picture 
of  these  easterly  waves  (see  Reed,  et  al,  1977  for  review). 

,It  is  known  that  most  tropical  cyclones  have  their  origins 
in  the  easterly  waves  (Gray,  1979).  Our  current  study  is  con¬ 
cerned  with  a  numerical  simulation  of  the  transformation  of 
easterly  waves  into  incipient  tropical  cyclones. 

4.2  THE  EASTERLY  WAVE  MODEL 

The  NRL/SAI  mesoscale  numerical  model  was  modified  to 
facilitate  the  simulation  of  easterly  waves.  For  this  simulation 
there  were  51  x  31  horizontal  grid  points  with  one  degree  re¬ 
solution  and  five  sigma  layers  in  the  vertical.  The  east-west 
lateral  boundaries  are  cyclic  and  the  second  derivatives  of 
thermodynamic  variables  normal  to  the  boundary  vanish  at  the 
north-south  boundaries.  Kuo's  (1975)  cumulus  parameterization 
scheme  and  Chang's  (1981)  boundary  layer  parameterization  scheme 
were  employed  in  the  model. 
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In  order  to  simplify  our  initial  attempt,  the  initial 
basic  easterly  current  is  set  uniformly  at  ^5  m  s'1.  The  Coriolis 
parameter  is  constant  throughout  the  model  domain  and  equals 

5  x  10" 5  s'1.  Observational  and  theoretical  studies  on  easterly 
waves  have  suggested  that  barotropic  (mainly  latitudinal  wind 

shear)  instability  and  beta  effects  are  important  for  the  intensification 
of  tropical  disturbances.  We  temporarily  leave  these  factors 
out  for  simplicity. 

4.3  PRELIMINARY  RESULTS 

The  initial  perturbation  is  introduced  by  a  prescribed 

6  h  stationary  heating  with  similar  spatial  distribution  as 
that  in  Section  2  but  with  a  lower  maximum  heating  rate  of 
20  K  day  ^  Fig.  4.1  shows  the  model  eastwest  cross-section 
of  latitudinal  components  of  the  wind  velocity  (v) ,  vertical 
velocity  (w),  and  relative  humidity  (RH)  at  48  h.  In  the  v 
field,  the  maximum  perturbations  occur  at  the  lowest  layer 

(o  »  0.975)  and  o  *  0.2.  The  amplitude  is  ^  n  s'1  with  cyclonic 
vorticity  it  low  levels  and  anticyclonic  vorticity  at  higher 
levels.  Compared  with  a  phase  III  GATE  composite  (Fig.  4.2, 
from  Reed  et  al ,  1977),  the  model  simulated  latitudinal  wind 
components  are  very  realistic  in  the  vertical  structure  and 
the  amplitude,  including  the  westward  tilt  of  the  zero  dividing 
line  of  south/north  components.  The  observed  maximum  perturbation 
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occurs  at  ^700  mb  instead  of  at  low  levels,  while  the  maximum 
vorticity  is  in  the  lowest  level  in  our  model.  This  is  related 
to  the  undesirably  strong  convective  heating  between  700  mb 
and  the  surface  in  our  model. 

The  simulated  vertical  motion  field  is  understandably 
noisier  than  the  observed  composite.  In  composite  analyses, 
small  scale  perturbations  such  as  gravity  waves  caused  by  moist 
convection  shown  here  in  our  model,  are  largely  smoothed  out. 
However,  the  maximum  upward  velocity  in  our  model  is  ^24  mb  h-*, 
four  times  stronger  than  the  observed.  We  recall  that  there 
are  no  barotropic  nor  barocllnic  instabilities  in  our  basic 
flow  field.  The  model  tropical  disturbance  is  mainly  sustained 
by  the  convection,  which  tends  to  produce  strong  vertical  motions 
in  spite  of  the  rather  realistic  v  components.  The  model  results 
can  be  improved  by  Including  barotropic  instability  in  the  basic 
easterly  flows. 

The  observed  location  of  the  maximum  upward  motion  is  about 
1/8- 1/4  wavelength  ahead  of  the  surface  trough  line,  whereas 
in  the  model  results  the  displacement  is  smaller.  This  decrepancy 
is  conjectured  to  be  the  result  of  the  absence  of  the  beta  effect 
in  the  model. 
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The  model  RH  field  agrees  well  with  the  observations.  There 
is  a  deep  moistened  troposphere  behind  the  trough  with  RH>70%  up 
to  ^400  mb.  The  midtroposphere  in  front  of  the  trough  is  slightly 
dryer  than  the  environment. 


4.4  SUMMARY 


Preliminary  results  of  an  attempt  to  simulate  easterly 
waves  are  presented.  In  this  attempt,  no  initial  barotropic 
nor  barfaclinic  instability  in  the  basic  flow  are  included. 

The  results  show  that  a  tropical  disturbance  is  realistically 
simulated  with  good  structure  in  the  horizontal  wind  field  and 
the  RH  field.  The  results  are  very  encouraging  and  suggest 
that  the  beta-effect  and  latitudinally  sheared  basic  easterlies 
be  included  in  future  attempts. 
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Fig  2.1  The  domain  of  the  monsoon  depression  model.  Contours 

are  surface  elevations  of  500  (dashed),  1000,  2000,  3000 


Fig  3.1 


The  horizontal  and  vertical  distribution  of  the  prescribed 
heating  pattern  used  in  the  Fujivhara  effect  experiments. 
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Fig  3.2  The  radial  distribution  of  the  quasi-steady  wind  speeds  gene¬ 
rated  by  the  stationary  prescribed  heating. 
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Fig  3. 


4  The  trajectories  of  stonn  centers  in  Exp.  1.  Numbers 
on  the  curve  denote  times  in  hour. 
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Fig  3.10  xhe  trajectories  of  the  free  drifting  storm  in  Exp.  7  (curve  C) 
and  of  two  Interacting  storms  in  Exp.  U  (curves  A  and  E). 


Fig  3.11  As  in  Fig  3.10,  except  for  Exps.  5  (solid  lines  with  dots) 
and  6  (solid  lines  with  squares). 
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Fig  4.1  East(right)-west(left)  cross-section  of  Che  easterly  wave 
model  showing  (a)  che  latitudinal  component  of  wind  vel¬ 
ocities  in  m  s”1  (upper  panel)  (b)  the  vertical  velocities 
in  w  h-1  (middle  panel),  and  (c)  the  relative  humidities 
in  X  (lower  panel)  associated  with  a  model  easterly  wave 
at  48  h. 


